Gauge Fields, Fermions and 
Mass Gaps in 6D Brane Worlds 



S. L. Parameswaran al , S. Randjbar-Daemi 62 and A. Salvio a3 



a International School for Advanced Studies, 
Via Beirut 2-4, 34014 Trieste, Italy 

b International Center for Theoretical Physics, 
Strada Costiera 11, 340 14 Trieste, Italy 



Abstract 

We study fluctuations about axisymmetric warped brane solutions in 6D minimal 
gauged supergravity. Much of our analysis is general and could be applied to other sce- 
narios. We focus on bulk sectors that could give rise to Standard Model gauge fields and 
charged matter. We reduce the dynamics to Schrocdinger type equations plus physical 
boundary conditions, and obtain exact solutions for the Kaluza-Klein wave functions 
and discrete mass spectra. The power-law warping, as opposed to exponential in 5D, 
means that zero mode wave functions can be peaked on negative tension branes, but 
only at the price of localizing the whole Kaluza-Klein tower there. However, remark- 
ably, the codimension two defects allow the Kaluza-Klein mass gap to remain finite even 
in the infinite volume limit. In principle, not only gravity, but Standard Model fields 
could 'feel' the extent of large extra dimensions, and still be described by an effective 
4D theory. 
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1 Introduction 



Minimal 6D gauged supergravity and its solutions have received much interest recently for 
several reasons. From the top down, the theory shares many features in common with 10D 
supergravity, whilst remaining relatively simple, and so it can be used as a toy model for 
10D string theory compactifications From the bottom up, it provides a context in which 
to extend the well-trodden path of 5D brane world models to codimension two. Moreover, in 
6D models, Supersymmetric Large Extra Dimensions (SLED) have shown some promise in 
addressing the two overshadowing fine-tuning problems of fundamental physics: the Gauge 
Hierarchy and the Cosmological Constant Problems 0. 

In terms of the phenomenological study of brane worlds, one should ask what are the 
qualitative differences between 5D and 6D models. For example, in five dimensional Randall- 
Sundrum (RS) models [SJ the warping of 4D spacetime slices is exponentially dependent on 
the proper radius of the extra dimension, whereas in the six dimensional models developed 
in Refs. [UIEI lEj it is only power law dependent. Moreover, the singularities sourced by the 
branes are distinct, the codimension one case being a jump and the codimension two case 
being conical. 

5D RS models with large (or infinite) extra dimension and the Standard Model (SM) 
confined to the brane were developed to explain the hierarchy in the Planck and Electroweak 
scales. Although the mass gap in the Kaluza-Klein spectrum goes to zero as usual in the 
infinite volume limit, 4D physics is retrieved thanks to the warp factor's localization of 
the zero mode graviton - and exponential suppression of higher modes - close to the brane 
with positive tension. Subsequently it was found that the step singularities in the geometry 
could also localize bulk fermions [7;, in much the same way as previously achieved with 
scalar fields and kink topological defects 8 . Models were then developed in which SM 
fields all 4 originate from the bulk as localized degrees of freedom [J2|. 

A study of warped brane worlds in 6D supergravity was intitiated in |H 151101) wn ere the 
focus was on the background solutions. As a result of supersymmetry ^H], and anomaly 
cancellation the 6D theory is highly constrained, and moreover fairly general classes 
of solutions have been obtained. In a general solution with 4D Poincare and 2D axial 
symmetry was derived, and it was shown that warping always leads to conical (or worse) 
singularities in the internal manifold, which can be naturally interpreted as 3-branes sources. 
It is certainly interesting to go beyond these background solutions, and study the dynamics 
of their fluctuations. The final objective would be to obtain the effective theory describing 
4D physics, and an understanding of when this effective theory is valid. 

Although in these constructions the SM is usually put in by hand, envisioned on a 
3-brane source, the bulk theory is potentially rich enough to contain the SM gauge and 
matter fields. As hinted above, should the SM arise as Kaluza-Klein zero modes of bulk 
fields, there are two ways to hide the heavy modes and recover 4D physics. They may have 
a large mass gap, and thus be unattainable at the energy scales thus- far encountered in our 
observed universe. Or they may be light but very weakly coupled to the massless modes, 
for example if the massless modes are peaked near to the brane, and the massive modes are 

4 Although the Higgs field should be confined to the brane in order not to lose the gauge hierarchy. 
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not. In any case, whether or not one expects the bulk to give rise to the SM, one should 
study its degrees of freedom and determine under which conditions they are observable or 
out of sight. This could also prove useful for a deeper understanding of the self-tuning 
mechanism of SLED, and its quantum corrections. 

A complete study of the linear perturbations is a very complicated problem, involving 
questions of gauge-fixings and a highly coupled system of dynamical equations. Some partial 
results have been obtained for the scalar perturbations in ^2j. In this paper we consider 
sectors within which SM gauge and charged matter fields might be found. By some fortune, 
these also happen to be two of the least complicated ones. Much of our discussion is general, 
and could easily be applied or extended to other 6D models with axial symmetry. We follow 
the usual Kaluza-Klein procedure, and reduce the equations of motion to an equivalent non- 
relativistic quantum mechanics problem, which we are able to solve exactly. We consider 
carefully the boundary conditions that the physical modes must satisfy, and from these 
derive the wave function profiles and complete discrete mass spectra. 

Our exact solutions enable us to analyze in detail the effects of the power-law warping 
and conical defects that arise in 6D brane worlds. We find that the warping cannot give rise 
to zero modes peaked at the brane, without also leading to peaked profiles for the entire 
Kaluza-Klein tower. On the other hand, the conical defects do break another standard 
lore of the classical Kaluza-Klein theory. Remarkably, even if the volume of the internal 
manifold goes to infinity, the mass gap does not necessarily go to zero. This decoupling 
between the mass gap and volume means that in principle SM fields, in addition to gravity, 
could 'feel' the extent of large extra dimensions, whilst still being accurately described by 
a 4D effective field theory. 

The paper is organised as follows. We begin in Section |21 by reviewing the 6D super- 
gravity theory, and its axially symmetric warped brane world solutions. Then, in Section 

01 we analyze the gauge field fluctuations, deriving the wave functions and masses of the 
Kaluza-Klein spectrum. A similar analysis is presented in Section 0] for the fermions. Sec- 
tion 13 discusses the physical implications of the results found, and in particular whether 
they can be naturally applied to the Large Extra Dimension scenario. Finally we end in 
Section El with some conclusions and future directions. 

In the appendices we give some results that are useful for the detailed calculations. Ap- 
pendix El sets our conventions, and Appendix iBl explains how the conical defects manifest 
themselves in the metric ansatz. In Appendix [0 we show in detail how the boundary con- 
ditions are applied to obtain a discrete mass spectrum, and we give the complete fermionic 
mass spectrum thus derived in Appendix iDl 

2 6D Supergravity and Axisymmetric Solutions 

Let us begin by presenting 6D J\f = 1 gauged supergravity, and its warped brane world 
solutions, whose fluctuations we will then study. The field content consists of a supergravity- 
tensor multiplet (Gmn, Bmn, cr, ipM, x) - with metric, anti-symmetric Kalb-Ramond field 
(with field strength Hmnp), dilaton, gravitino and dilatino - coupled to a number of gauge 
multiplets (Am, A) - with gauge potentials and gauginos - and a number of hypermultiplets 
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(<£, ^) - with hyperscalars and hyperinos. The fermions are all Weyl spinors, satisfying 
^TipM = V'M) ^7X = ~Xi rVA = A and = —*ff, and in general the theory has anomalies. 
However, for certain gauge groups and hypermultiplet representations these anomalies can 
be cancelled via a Green-Schwarz mechanism [111 I13| . We will consider a general matter 
content, with gauge group of the form Q = Q x U(1)r, where U(1)r is an R-symmetry, and 
Q is in general a product of simple groups. For example we could take the anomaly free 
group Q = Eq x Ej x U(1)r, under which the fermions are charged as follows: ipM ~ (1, 1) 1 , 
X ~ (1, l)i, A ~ (78, 1)! + (1, 133)! + (1, l)i, * ~ (1, 912) dJ. 
The bosonic action takes the form 5 |1U| 



S B 
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where k represents the 6D Planck scale, g and g\ are the Q and U(l)p gauge coupling 
constants respectively 6 and Fmn arid F\mn the corresponding field strengths. G a p($) is 
the metric on the target manifold of the hyperscalars, and here the index a runs over all 
the hyperscalars. The dependence of the scalar potential on $ Q is such that its minimum 
is at <1> Q = 0, where it takes a positive-definite value, v(0) = gf, due to the R-symmetry 
gauging [141 115j . We therefore fix <£ a = 0. The remaining equations of motion (EOM) are 
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where Dm is the gauge and Lorentz covariant derivative. 

We will consider a general class of warped solutions with 4D Poincare symmetry, and 
axial symmetry in the transverse dimensions: 



ds< 



G MNdX^ dX^ 
A 
a 

Hmnp 



e^^dx^dx" + dp 2 + e^ p >d(p 
A v (p)Qd(p, 

0, 



B(P), 



(2.3) 



5 We give conventions and notation in Appendix^ For fermionic terms see |10|. 

6 In general, since Q consists of several simple factors, g represents a collection of independent gauge 
couplings. 
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with < p < ~p and < (p < 2ir. Here p, v = 0,1,2,3, A is a gauge field, and Q is a 
generator of a U(l) subgroup of a simple factor of Q, satisfying Tr (Q 2 ) = 1. 
In the following we shall also use the radial coordinate defined by 

u(p) = [ P dp'e- A ^ 2 , (2.4) 
Jo 

whose range is < u < u. In this frame the metric reads 

ds 2 = e A ^ (r jflu dx^dx u + du 2 ^j + c B{ ^d^ 2 . (2.5) 



Given the above ansatz, the general solution to the equations of motion (|2.2|) has been 
found by Gibbons, Giiven and Pope (GGP) in [I] 7 . Although much of our formalism for the 
perturbation analysis can be applied to the general ansatz (J2.3|) . we will focus on a subset 
of this general solution, namely that which contains singularities no worse than conical. 
Thus, in addition to the ansatz H2.3|) . we impose the following asymptotic behaviour for the 
metric: 



e A p —>° constant / 0, e A p —> p constant ^ 0, 



and 

e BP-_ 



(l-5/27r) 2 p 2 , e BP ^(l-5/27r) 2 (p-p) 2 , (2.6) 

that is we assume conical defects with deficit angle S at p = and <5 at p = p, at which 
points the Ricci scalar contains delta-functions (see Appendix lB|) . These can be interpreted 
as 3-brane sources with tensions T = 25/ 'k 2 and T = 25/ 'k 2 |16j . 
The explicit conical-GGP solution is then 8 [3]: 

e A = e na/2 = [E ^ _ j COt 2 (u/r ) 

V fo ' fl 

A = -^rQdtp, (2.7) 

where a is a real number, which we can take to be positive without loss of generality, q is a 
real number and g is the gauge coupling constant corresponding to the background gauge 
field. For example, if A lies in U(1)r, then g = g±. Also, 



/o-l + cot 2 (-), /i-l + 5^ 2 (-), (2. 



„2 
_0 

So/ r\ \r 

with rg = k 2 /(2g 2 ), r\ = 8/q 2 . Moreover, u = irro/2 and the deficit angles are 



5 = 2vr ^1 - a-ij , (2.9) 
5 = 2tt(1 -a) . (2.10) 



7 Relaxing the condition of 4D Poincare symmetry to that of only 4D maximal symmetry should allow 
more general solutions to be found |17) . More general solutions also exist breaking axial symmetry or 
having nontrivial VEVs for the hyperscalars |19| . 

8 The coordinate u is related to the coordinate r in 4 by r = ro cot(it/ro). 
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Notice then that a non-trivial warping enforces the presence of a 3-brane source. On the 
other hand, the parameter a is not fixed by the EOM and it represents a modulus. 

The expression for the gauge field background in eq. (|2.7|) is well-defined in the limit 
u — > 0, but not as u — > u. We should therefore use a different patch to describe the 
u = u brane, and this must be related to the patch including the u = brane by a single- 
valued gauge transformation. This leads to a Dirac quantization condition, which for a field 
interacting with A through a charge e gives 

4< *0 n 9 AT to 1 1 \ 

e = ea = A , (2.11) 

nq r gi 

where N is an integer. The charge e can be computed once we have selected the background 
gauge group, since it is an eigenvalue of the generator Q. Finally, note that the internal 
manifold corresponding to solution (|2.7|) has an S 2 topology (its Euler number equals 2). 

We end this section by considering the various parameters in the model, and the phe- 
nomenological constraints which can arise when we give it a brane world interpretation. 
There are three free parameters in the 6D theory, which can be taken to be the gauge cou- 
pling g, and two out of the following three parameters: the 6D Planck scale, k, the gauge 
coupling g\ and the length-scale tq = K,/v2g\. In the solution there are two free parameters, 
r% (or q) and a. However, one combination of all these parameters is constrained by the 
quantization condition (|2.11|) . 

The relation between the 6D Planck scale n and our observed 4D Planck scale K4 is 

^ V2 = h (2 - 12) 

where the volume V2 is given by 

V 2 =[ d 2 yV^Ge~ A = 2tt / due^ A+B ^ 2 . (2.13) 
For solution (|2.7|) we have 

V 2 = 4ma{^\ 2 . (2.14) 

Notice that this volume does not depend on n, and so we can keep it fixed whilst varying 
the warp factor, namely e A in (|2.7|) . Moreover from (|2,12|) a phenomenological constraint 
follows between the bulk couplings and the brane tensions, which can be written: 

y= = J-K*. (2.15) 



This implies that g\j \fa is very small, of the order of the Planck length. 

Now let us embed the Arkani Hamed-Dimopoulos-Dvali (ADD) scenario [20] into the 
present model, in order to try to explain the large hierarchy between the Electroweak 
scale and the Planck scale via the size of the extra dimensions. Thus identifying the 6D 
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fundamental scale with the Electroweak scale n ~ TeV 2 and constraining the observed 4D 
Planck scale n\ ~ 10~ 30 TeV~ 2 ; the above relation translates to 9 : 

y/aro ~ 0.1mm. (2.16) 

Here, the large extra dimensions corresponds to tuning the bulk gauge coupling and brane 
tensions. However, we can also observe that (j2.16|) fixes just one parameter among a, ro and 
r± and we still have two independent parameters even if we require large extra dimensions. 
Later we will see that this novel feature proves to have interesting consequences for the 
mass spectrum of fluctuations. 

3 Gauge Fields 

Having established the brane world solution and its properties, we are now ready to examine 
the fluctuations about this background, which will represent the physical fields in our model. 
In this section our focus will be on the gauge field fluctuations. 

Normalizable gauge field zero modes in axially symmetric codimension two branes are 
known to exist |21 |l22| l2*H]. However, in these known examples there is no mass gap between 
the zero and non-zero modes which renders an effective 4-dimensional description somewhat 
problematic, especially in non-Abelian case [21]. I n contrast to this for the axisymmetric 
solutions studied in this paper the presence of a mass gap will be automatic due to the 
compactness of the transverse space. In this section we shall give the full spectrum of zero 
and non-zero modes. 

Given the symmetries of the problem, we can expect that the gauge fields in the low 
energy effective theory belong toHx U(1)kk, where H is the unbroken subgroup of Q that 
commutes with the U(l) C Q in which the monopole lies 10 . The U(1)kk arises from the 
vector fluctuations of the metric, due to the axial symmetry of the internal manifold, and 
is promoted to SU (2)kk in the sphere limit of the background. 

The non-Abelian sector of Ti. may be rich enough to contain the SM gauge group. For 
example, consider the anomaly free model of [TT| . with gauge group Q = Eq x E-j x U(1)r, 
and the monopole background in Eq. The surviving gauge group, SO(W) x Ej x U(1)r x 
U{1)kKi then contains the Grand Unified Group SO(10), and the model also includes 
charged matter in the fundamental of 50(10). Therefore our present interest will be in the 
fields belonging to various representations of TL. Specifically, we will consider gauge field 
fluctuations orthogonal to the monopole background. For the case Q = Eq x Ej x U(1)r, 
the gauge field sectors that are covered by our analysis are given in Tabled 

3.1 Kaluza-Klein Modes 

Using the background solution in the 6D action (|2.1|) , we can identify the bilinear action for 
the fluctuations. This step requires some care, because to study the physical spectrum we 

We use the following conversion relation: (TeV) -1 ~ 1CT 16 mm. 
10 Due to the Chern-Simons coupling in supergravity, the (7(1) gauge field in the direction of the monopole 
eats the axion arising from the Kalb-Ramond field and acquires a mass |1 1 1 IT], 
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must first remove the gauge freedoms in the action due to 6D diffeomorphisms and gauge 
transformations. The problem has been studied in a general context in [23], where the 
authors choose a light-cone gauge fixing. 

In the light-cone gauge, the action for the gauge field fluctuations, orthogonal to the 
monopole background, at the bilinear level reads [23] 

s G (v, v) = -J d 6 xV^G ^ (dfjyjdrv* + e- A d p Vjd p v 3 + d v v 3 d^) , (3.1) 

where Vj is the gauge field fluctuation in the light cone gauge (j=l,2) and all the indices 
in 1)3.1 J) are raised and lowered with the /^-dependent metric Gmn given in ()2.3j) . Indeed, 
here and below Gmn represents the background metric. We have multiplied the formula of 
:23 by an overall e^, with <j> = kct/2 and a in the background, due to the presence of the 
dilaton in our theory 11 . Notice that since we are looking at the sector orthogonal to the 
monopole background, the Chern-Simons term does not contribute, and the action takes a 
simple form. 

In general, the covariant derivative D^Vj includes the gauge field background 

DtpVj = d v Vj + ieyApVj, (3.2) 

where again the charge ey can be computed using group theory once the gauge group Q is 
chosen. The value ey = corresponds to the gauge fields in the 4D low energy effective 
theory. However, since we can do so without much expense, we keep a generic value of ey. 
Those fluctuations with ey 7^ corresponds to vector fields in a non-trivial representation 
of the 4D effective theory gauge group. The Dirac quantization condition ()2.11)) then gives 
ey Aag/(Kq) = Ny, where Ny is an integer. 

Next we perform a Kaluza-Klein expansion of the 6D fields. Since our internal space is 
topologically S 2 , we require gauge fields to be periodic functions of ip: 

Vj(X) = ]T V jm (x)f m (p)e im *, (3.3) 

m 

where m is an integer. 

If we put ()3.3)) in 1)3.1)1 we obtain kinetic terms for the 4D effective fields proportional 

to 

f d A xY,r}^d„V] m d v V im f dpe* +B / 2 \f m \ 2 . (3.4) 

1 That the dilaton invokes only this simple change with respect to Ref. |2o| can be seen as follows. First, 
notice that since we are considering fluctuations orthogonal to the U(l) background, there are no mixings 
with other sectors, and the bilinear action is simply S G = -1/4 / <f X*/^Ge' t 'G MN G PQ (F mp F N q) {2) . We 
emphasise that Gmn and (f> now signify the background fields. Also, 0' 2 - 1 indicates the bilinear part in the 
fluctuations. Next, make the change of coordinates, dp — e -< ^ 2 d/5, and rewrite the background metric in 
as ds 2 = e -0 (e^rj^dx^dx" + dp 2 + e^dip 2 ^J , with A= A + <f> and B = B + 4>. In this way, the bilinear 

action, Sg, reduces to exactly the same form as that of Ref. |23j . and we can proceed as they do to transform 
into light-cone coordinates, fix the light-cone gauge, and eliminate redundant degrees of freedom using their 
equations of motion. 
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Therefore physical fluctuations, having a finite kinetic energy, must satisfy the following 
normalizability condition (NC): 

du\ip\ 2 < oo, (3.5) 



where 

V = e^ +A+B ^f m . (3.6) 

The quantity \ip\ 2 represents the probability density of finding a gauge field in [u,u + du]. 

In fact, this is not the only condition that the physical fields must satisfy. If we want to 
derive the EOM from (|3.1j) through an action principle we have to impose 12 the following 
boundary condition HH1 HZ] 

J d 6 Xd M (^Ge^ A Vp M V^ = , (3.7) 

where T>m is the gauge covariant derivative. Equation 1)3 . 7 j) represents conservation of 
current Jm = VjT>MVj. Moreover, since the fields are periodic functions of <p, (|3,7j) 
becomes 

" e <i> ~ A VjdpVj] P = 0. (3.8) 



The EOM can then be derived as: 

V^Ge^- 2A rj^d^d u Vj = -d p (V^Ge^dpVj) - v^Ge^^D 2 Vj. (3.9) 
By inserting ([3.3)1 in 1)3.9)) we obtain 

-</>+2A 

- -jT£ dp {^-^ e<t " Ad pfm) + e A ~ B (m + e v A v f f m = M^ m f m , (3.10) 

where My m are the eigenvalues of rj^dpdy. 

At this stage, we can already identify the massless fluctuation that is expected from 
symmetry arguments. For ey = 0, when m = 0, a constant /o is a solution of 1)3.10)1 with 
My = 0. This solution corresponds to 4D effective theory gauge fields. It has a finite 
kinetic energy, and trivially satisfies (|3.8|) . The fact that such gauge fields have a constant 
transverse profile guarantees charge universality of fermions in the 4D effective theory (see 
below) . 

To find the massive mode solutions, we can express 1)3.10)1 in terms of u and ip and 
obtain a Schroedinger equation: 

(-dl + V)^ = M 2 i;, (3.11) 

where the "potential" is 

V{u) = e A - B (m + e v A^) 2 + e -(2^+B)/4gg e (2*+yi + B)/4_ (3 12) 

12 Actually we impose that for every pair of fields Vj and V- the condition / d 6 Xd M ( % /^Ge 0_A V}£> M V-) = 
is satisfied but in 13.71 the prime is understood. 
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We want to find the complete set of solutions to 1)3.11(1 satisfying the NC 1)3.5)) and the 
boundary conditions ()3.8|) . which can be written in terms of u and tp as follows 

( lim - lim) U* \-d u + \ (2d u cf> + 8 U A + d u B)} A = 0. (3.13) 

In order for ((3.13)) to be satisfied, both the limits u — > and u — ► « must be finite. Condition 
1)3.13)1 ensures that the Hamiltonian in the Schroedinger equation 1)3.11)1 is hermitian, and 
so has real eigenvalues and an orthonormal set of eigenfunctions. Therefore, we shall call it 
the hermiticity condition (HC). 

So far our analysis has been valid for all axially symmetric solutions of the form 1)2.3)1 . 
We will now use these results to determine the fluctuation spectrum about the conical-GGP 
solution ()2.7j) . We observe that V(u) then contains a delta-function contribution, arising 
from the second-order derivative of the conical metric function d 2 B (see AppendixlBland eq. 
I)B.6)I ). However we can drop it because d 2 B also contains stronger singularities at u = 
and u = u: respectively 1/u 2 and l/(u — u) 2 . These singularities are also a consequence of 
the behaviour of e B given in 1)2.6)1 and they imply that the behaviour of the wave functions 
close to u = and u = u cannot depend on the mass. In particular, this immediately implies 
that if the wave functions of zero modes are peaked near to one of the branes, then the same 
will be true also for the infinite tower of non-zero modes. In other words, we cannot hope 
to dynamically generate a brane world scenario, in which zero modes are peaked on the 
brane, and massive modes are not, leading to weak coupling between the two sectors 13 . If 
we are to interpret the zero mode gauge fields as those of the SM, therefore, for the massive 
modes to have escaped detection they must have a large mass gap. 

Meanwhile, we note that in contrast to the non-relativistic quantum mechanics problem, 
here we cannot deduce qualitative results about the mass spectrum from the shape of 
the potential. This is because the boundary conditions to be applied in the context of 
dimensional reduction are in general different to those in problems of quantum mechanics. 
In particular, the HC ()3.13)) is a non-linear condition, contrary to the less general linear 
boundary conditions usually encountered in quantum mechanics to ensure hermiticity of 
the Hamiltonian. We will be able to impose the more general case thanks to the universal 
asymptotic behaviour of the Kaluza-Klein tower. 

Returning then to our explicit calculation of the Kaluza-Klein spectrum, we can write 
V(u) as 



u 



and 



V(u) = V + v cot 2 ( — J + vtan 2 f — 1 , (3.14) 

TqVq = 2raw(m — Ny)TJj — -, TqV = m 2 uj 2 — -, r^v = (m — Ny) 2 uJ 2 — —. (3.15) 
Moreover in this case the expression ()3.13)) for the HC becomes 

lim V/ (~d u + l^-z) ^ ~ lim V* (~d u + V = 0. (3.16) 
u >u \ 2 u — u J u— »Q \ 2u J 

13 In fact, a similar singular behaviour for the potential in general arises for the general axisymmetric 
solutions given in |lj and studied in [B], where the hypothesis 12.61 is relaxed. 
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If we introduce z and y in the following way |2H] 

i; = z^(l-zfy(z), (3.17) 

equation (|3.11|) becomes 

z(l - z)d\y + [c - (a + b + l)z] - a&y = 0, (3.18) 

where 

7 = i [1 + 2(m - JW)5?] , (3=^(1 + 2mw) , c = 1 + (m - N v )uJ, 

a = \ + ^ + V m ~ N v^ + \"J r o M v,m + 1 + [muj - (m - N v p] 2 , 

b = ~ + j u + ~(m - iV v )u7 - ~ m + l + [ m u-(m- N v )u} 2 , (3.19) 

and 

u = (l-5/2ir)- 1 , w = (1 - 8/2-n)- 1 , (3.20) 

Equation (|3.18|) is the hypergeometric equation and its solutions are known. For c 7^ 1 the 
general solution is a linear combination of the following functions: 

yi (z) = F(a,b,c,z), y 2 (z) = z 1 - c F(a + l-c,b + l-c,2-c,z), (3.21) 

where F is Gauss's hypergeometric function. So for c / 1 the general integral of the 
Schroedinger equation is 

i> = K x i>x + K 2 ifo, (3.22) 

where 

i/ji = z\l - zf yi . (3.23) 

and K\p. are integration constants. For c = 1 we have ip\ = ip 2 but we can construct a 
linearly independent solution using the Wronskian method and the general solution reads 

/« d u ' 
TITTY ( 3 - 24 ) 

Now we must impose the NC (|3.5|) and HC (j3.16|) . to select the physical modes. In 
Appendix Owe give explicit calculations; the final result is that the NC and HC give the 
following discrete spectrum. The wave functions are 

ip oc z 7 (l - zfF{a, b, c,z), for m>N v , (3.25) 
ip cx z 7+1 ~ c (l - zfF{a + 1 - c,b + 1 - c,2 - c,z), for m<N v . (3.26) 

and the squared masses are as follows: 
• For N v <m<0 

M Vn,m = ^ [n(n + 1) + Q + r») [-mu> + (m - N v p] } > 0. (3.27) 
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For m > Ny and m > 



M Vn,m = ^2 + l ) + Q + n ) [ mU; + ( m ~ Wy)E7] + mw(m - JVy)wJ > 0. 

(3.28) 



• For m < AV and m < 

4 



n(n + 1) + (— + n^j [—mco + (Ny — m)u] — muj(Ny — m)u\ > 0. 

(3.29) 



For < m < Ny 

M Vn,m = ^2 [ n ( n + !) + Q + n ) I mw + (^V " m)57]| > 0. (3.30) 



The masses given in (|3.27|) and ()3.28|) correspond to the wave function (|3.25j) whereas the 
masses given in (|3,29[) and ()3.3U|) correspond to the wave function ()3.26|) . We observe that 
there are no tachyons and that the only zero mode is for n = 0, m = and Ny = (ey = 0), 
corresponding to gauge fields in the 4D low energy effective theory. 

As a check, we can consider the S 2 limit (lo,uj — > 1), whose mass spectrum is well-known. 
Our spectrum (|3.27|) - ([3.30[) reduces to 

a 2 My = 1(1 + 1) - (~y) i multiplicity = 21 + 1, (3.31) 

where 14 I = \^-\ + k and k = 0, 1, 2, 3, .... This is exactly the result that one finds by using 
the spherical harmonic expansion [22] from the beginning. 

At this stage we can point towards a novel property of the final mass spectrum. Observe 
that in the large a (small ZJ) limit the volume V2 given in equation (|2.14|) becomes large 
but the mass gap between two consecutive Kaluza-Klein states does not reduce to zero as 
in standard Kaluza-Klein theories 15 . This a consequence of the shape of our background 
manifold and in particular of the conical defects. Notice that the large a limit corresponds 
to a negative tension brane at u = u, but not necessarily at u = 0. 

In Section 0] we will show that the same effect appears also in the fermionic sector, and 
we will turn to a discussion of its implications in Section |3 



3.2 4D Effective Gauge Coupling 

Let us end the discussion on gauge fields by briefly presenting the 4D effective gauge cou- 
pling. This can be obtained by dimensionally reducing the 6D gauge kinetic term. We 

14 The number I is denned in different ways in equations I3.27p - H3.30p . For instance we have / = n+\Nv /2| 
for $FF7\ . 

15 This is also true for the proper volume of the 2D internal manifold. 
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consider the zero mode fluctuations in Tt, about the background (|2.3|) . so that the 4D 
effective gauge kinetic term is: 



d b Xy 



d ko/2 



TrF MN F 
1 



MN 



dud^ A+B V 2 ft 



Recalling that = const and normalizing it to one, we can read: 

1 V 



9eff 



(3.32) 



(3.33) 



4 Fermions 

We will now consider fermionic perturbations, and in particular our interest will be in 
the sector charged under the 4D effective gauge group, Tt, discussed above. These fields 
arise from the hyperinos and the Tt gauginos, for which we also restrict ourselves to those 
orthogonal to the U(1)r- Thus we are considering matter charged under the non-Abelian 
gauge symmetries of the 4D effective theory. For instance, for the anomaly free model 
Eq x E-j x U(1)r, with the monopole embedded in the Eq, the gauginos in the 78 of Eq 
contain a 16 + 16 fundamental representation of the grand unified gauge group 50(10), 
and our analysis will be applicable to them. In TabelQ we & ve the complete list of fermion 
fields that are included in our study, for the said example. 

We proceed in much the same way as for the gauge field sector of the previous sec- 
tion, transforming the dynamical equations and necessary boundary conditions into a 
Schroedinger-like problem, to obtain the physical modes and discrete mass spectrum. 

The bilinear action for the fluctuations of interest takes a particularly simple form, 
comprising as it does of the standard Dirac action: 



5 p 



d b Xy 



-G \T M D M \ 



where 16 



D M \ = ( d M + ±-Q [ m' B] [T A , T B ] + ieA M ) A 



(4.1) 



(4.2) 



Here e is the charge of A under the U{\) monopole, and Gmn, ^m'^' anc ^ Am are the back- 
ground metric, spin connection and gauge field corresponding to an axisymmetric solution 
(|2.3j) . Analogously to the gauge field analysis, in order to derive the Dirac equation 

F M D M X = (4.3) 

from (|4.1|) by using an action principle (3111) we require conservation of fermionic current: 

/ d 6 Xd M (V^G XT M \) = 0. (4.4) 



6 Our conventions for T A and fi^' S ' are given in Appendix 1X1 
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Equation 1)4.4)1 implies that the Dirac operator F M Dm is hermitian, and we shall again refer 
to it as the HC. Our aim is to find the complete fermionic spectrum, that is a complete set 
of normalizable solutions of 1)4.3(1 satisfying 1)4.4(1 . 

Some care is now needed when discussing the background felt by the fermionic sector 
in (|4.2|) . As already mentioned, in order to have a correctly defined gauge connection, it is 
necessary to use two patches related by a single- valued gauge transformation. The same is 
true for the spin connection, which must be defined in such a way as to imply the conical 
defects in the geometry. Henceforth we focus on the patch including the p = brane, chosen 
to be < p < p. For this patch a good choice for the vielbein is 



jx _ e A/2,a < a \ _ ( cos V -e B/2 shnp 

S-e V \e m )- I gin ^ e B/2 CQS(p 



(4.5) 



where a is a 4D flat index, a = 5, 6 a 2D flat index and m = p, if. The corresponding spin 
connection is 

1 

~ /! ~ 2" 

n = Ol 5 ' 61 = (l- \B'e B / 2 ) , (4.6) 



0. 



where ' = d p . It can be checked that this gauge choice correctly reproduces Stokes' theorem 
for a small domain including the conical defect . 

We are now ready to study the Dirac equation 1)4.3(1 for 6D fluctuations, and write it in 
terms of 4D effective fields. Since A is a 6D Weyl spinor we can represent it by 



A 



Ai 




(4.7) 



where A4 is a 4D Dirac spinor: A4 = Xr + Xl, 7 5 Aj? = Xr, 7 5 Al = — Xl- By using the ansatz 
1)2.3)) . the vielbein 1)4.5)1 . the spin connection 1)4.6)1 and our conventions for T A in Appendix 
lAl the Dirac equation 1)4.3(1 becomes 



e- A / 2 ^d,X R 



vp 



-dp — ie 



-B/2 



{d v + ieA v ) -A' + -ne~ B / 2 



-vp 



d p - ie- B ' 2 {8 V + ieA^) + A' - \sie' B / 2 



Xr, 
Xl- 



Performing the Fourier mode decomposition: 

A 4 (A) = Xr(X) + X L (X) = £ (X R , m (x)fR, m (p) + X L , m (x)f L>m (p)) e im ^, 



(4.8) 
(4.9) 

(4.10) 



where m is an integer, and inserting into 1)4.8)1 and 1)4.9)1 we find: 



e A/2 ^dpX L:m+1 f L , m+ i 
e~ A/2 ^dpXR, m -ifR, m -i 



dn + e 3/2 { m + -SI + eA u 



o + (T BI2 (m - io + eA^] + A' 



-4 Aij )m /i? jm ,(4.11) 

XL,mfL,m- (4.12) 



7 See Appendix FBI for some steps in this calculation. 
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For the boundary conditions, analogously to the gauge fields, the NC can be found to 



be: 



where 



du < oo 



(4.13) 



* = e A+B ^f Rm (4.14) 
and a similar condition for left-handed spinors. Meanwhile, the HC (|4.4j) can be written: 



R,m 



0. 



(4.15) 



Having set up the dynamical equations and the relevant boundary conditions, we shall 
now use this information to study the complete fermionic spectrum in Subsections 14.11 and 
14.21 In particular, we will focus on the questions of wave function localization, and the mass 
gap problem, crucial to the development of a phenomenological brane world model. 



4.1 Zero Modes 

We begin by finding the zero mode solutions, for which the problem simplifies considerably. 
Indeed, for the zero modes 7^<9^ = 0, and the equations for right- and left-handed modes 
(|4~TT|) and (|4~T2l decouple: 



d p -e- B/2 (m 



1 



+ eA v )+A' --n e - B ^ 2 



fR,m = 0, 
fL,m = 0. 

By using the expression for Q in equation Q4.6JI . the solution of Q4.16JI is 



d p + e- B ' 2 {m + eAp) + A' - \^e~ B / 2 



fR,m(p) oc exp 



' A ~\ B + f d p' e ' B/2 ( m + \ + e A 



(4.16) 
(4.17) 

(4.18) 



whereas the solution of Q4.17JI can be obtained by replacing m,e — * —m,—e in Q4.18JI . 
The solution ()4.18j) for e = was found in [SJ- Here we give the expression for every e 
because we want to include charged fermions. We note that the zero mode solution (|4.18f) 
automatically satisfies the HC given in (|4.15jl . From (j4,18|) . (|4.14|) and (|2.4|) we obtain 



ip oc exp 



<!«'<< [ - U) 2 {m + ^+eA^ 



For the conical-GGP background, (|2.7|) . the explicit expression for ifi is 



u s w(l/2+m) / u X S7(JV-l/2-m) 

ip oc sin ( — I cos — 

To/ Vro. 



where we used (|2.11(l . The NC is satisfied when 



5 5 
1 < m < iV - — . 

4-7T 47T 



(4.19) 



(4.20) 



(4.21) 
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From here we retrieve the result that for the sphere, which has 5 = 5 = 0, there exist 
normalizable zero modes only for e ^ (N 7^ 0), that is for a non- vanishing monopole 
background [22]. Moreover, as found in [HJ, we see that the conical defects also make 
massless modes possible, provided that there is at least one negative deficit angle, even if 
N = 0. However, (|4.21|) implies that for positive tension branes, 5, 5 > 0, the adjoint of 
Tt, which has e = 0, is projected out. If 7i contains the SM gauge group, this is appealing 
since the fermions of the SM are not in adjoint representations. In any case, the number of 
families depends on 5, 5 and N. 

Let us now consider the wave function profiles, (|4.2()jl . Observe that ip is peaked on the 
u = brane (that is, ip — ► oo as u — > 0, and ip — > as u — > u) when 

m<-l/2, and m<-l/2 + JV. (4.22) 

By comparing (|4.21j) and Q4.22JI we understand that we have normalizable and peaked ip 
only for 5 < (that is for negative tension brane). If N — 5/4tt > we can have normalizable 
zero modes for 5 > (positive tension brane) but the corresponding ip are not peaked on 
the u = brane (indeed, ip — * as u — > 0). On the other hand from (|4.2U|) we see that ip is 
peaked on the u = u brane when 

m>-l/2, and m> -1/2 + N. (4.23) 

By comparing (|4.21j) and (|4.23j) we understand that we have normalizable and peaked ip 
only for 5 < (that is for negative tension brane). 

We can also analyze the chirality structure. Since the left handed wave functions can 
be obtained from (j4,20|) by replacing m, N — > —m, —N, in order for them to be peaked on 
the u = brane we need 

m>l/2, and m> 1/2 + N, (4.24) 

whereas in order for the left handed wave functions to be peaked on the u = u brane we 
need 

m<l/2, and m< 1/2 + N. (4.25) 

So, if we were to require that ip be peaked on a brane, we always have a chiral massless 
spectrum because the chirality index counts the difference of modes in f^ m and fh,m with 
given m, Nn{m) — NL{m) |31j . We should point out, however, that in fact peaked zero 
modes, ip, may not be necessary in order to have an acceptable phenomenology. The answer 
to this question can be found only after constructing the complete spectrum, and studying 
the couplings between different 4D effective fields. 

4.2 Massive Modes 

We now move on to a study of the complete Kaluza-Klein tower for the fermions. We begin 
by establishing the corresponding Schroedinger problem. The two coupled first order ODEs, 
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equations (|4.11j) and 1)4. 12|) . can be equivalently expressed as a single second order ODE 
and a constraint equation, as follows 18 



-d p + hd p + g m ) f R>rn = M Fm f R 

MF,mfL,m+l = e A ^ 2 

where M Fm are the eigenvalues of {^d^) 2 and 



d p - A 1 + (m + ~0 + eAj) e~ B ' 2 



fR,', 



h = -^A' + (n-l)e~ B / 2 , 



In/ 1 



m 



5 A „ 



m 



-Vl' - -VLB' - —B' + -A'VL + — - 1 )A' - -B'A^ + eA' + -A' A 



(4.26) 
(4.27) 



(4.28) 

-B/2 



+ 



1 VL 2 

m(m + 1) + -VL - — + (2m + l)eA^> + e^; 



2 42 



-B All ^ / 4/\2 



< - - A" - 2 . (4.29) 



Once /jt, m is known we can compute fL,m+l by using Q4.27JI . so we can focus on and 
study the second order ODE (|4.26|) . If we express this equation in terms of if) and u we 
obtain the Schroedinger equation 



(-d 2 u + v)^ = M 2 ^, 
where the "potential" V is given by 

V(u) 



(4.30) 



ed u A^ A - B ^ 2 + (~ + m + eA^j d u e^l 2 



+ 



1 2 

- + m + eA<p + (m + eA^) 



(4.31) 



We observe that transformation ()4.14j) exactly removes the delta-functions which appear in 
(|4.29|) through Vt' (see Appendix [Bj). However, just as for the gauge fields, a singular be- 
haviour is observed in the potential, so that the asymptotic behaviour of the wave functions 
does not depend on their mass. 

Our problem is now reduced to solving equation ()4.3U|) with the conditions NC (|4.13f) 
and HC (|4.1S|) . By using (|4.27j) for M F ^ and definitions (|2~i|) and (|4.14j) we can rewrite 
(14.151) as follows 



lirn - lim) V* \-d u + (m + \ + eA v ) e (A ^ )/2 



u^u u—*0/ 



ip = 0. 



(4.32) 



We can now proceed in exactly the same way as for the gauge field sector. For the 
conical-GGP solution ()2.7|) the explicit expression for V has the form ()3.14j) . but now with 



r 2 V 



l - + m 



jj — x + 2^'^ - [ — + m — N 



(4.33) 



18 Ea, 14.26^ can also be obtained by squaring the 6D Dirac operator, and using that [Dm,Dn] = 



\R mn ab Yab +ieFMN- 
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+ to 
+ to 



— U! + U) 
ZJ + UJ 2 



+ TO 



+ m - 2N 



+ UN (uJN-1). 



Moreover, in this case the explicit expression for the HC is 



lim ifj* 



n _m + l/2-N 
-O u + U — 



u — u 



ib — lim ib* 



-d u + LO 



m + 1/2 



ip = 0. 



As in the gauge fields sector we introduce z and y in the following way 

ij) = (1 - z f y(z), 



z = cos I — 

To 



(4.34) 
(4.35) 

(4.36) 

(4.37) 



so that equation (|3.14f) becomes a hypergeometric equation (|3.18|) . with parameters: 



7 



1 + tJ I - + TO- iV 



1 

2 _ 

lJ /l 
1 + - - + TO- iV 

2 \2 
uJ /l 



/3= t 1-w 



+ TO 



+ tJ I - + TO - iV 



1 + — [ - — //) - .V ] - — ( - + m 

„2 7\ #2 



2 V2 



2 \2 



A = rgM£ m + (UN) 2 + ( | + m 



1 



1 



lJ(uJ - 2u)[-+m — N 



- uj 2 N + u 2 ( - + TO 



We can construct two independent solutions and ^ of the Schroedinger equation (|3.14|) 
as in Section |S1 and impose the NC (|4,13|) and the HC (|4.3fi[) to obtain the physical modes. 
The resulting wave functions are: 



(4.38) 



where the integration constants, i^i,2, are fixed in Appendix[D] We plot a few of the wave 
function profiles in Figure ^ 

The complete discrete mass spectrum is also given in Appendix[DJ There it can be seen 
that the same finiteness of the mass gap in the large a (hence large volume) limit, found in 
the gauge field spectrum, can be observed here. Moreover, for a ~ 1, the mass gap between 
the zero modes and the massive states now goes as: 



M, 



GAP 



r 2 + r 2 



(4.39) 



Therefore, for the fermions, a finite mass gap in the large volume limit can also be obtained 
by taking tq —* oo and turning on 8, thus allowing r% to remain finite. This contrasting 
behaviour to the standard Kaluza-Klein picture is a consequence of the conical defects 
(uJ 7^ 1 and to ^ 1) in our internal manifold. Below we shall consider its implications for 
phenomenology. 
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Figure 1: Fermion Wave Function Profiles: n — 0,1,2 modes plotted for angular momentum numbers 
rn — —1,0 (eqs l|D.6(l and respectively). The parameters are chosen to be (ro,u,Z3, e) = (1, 1/4, 1,0), 

corresponding to a single negative tension brane at u = 0. Also the normalisation constant is set to 1. 
The number of intersections with the it-axis equals n, according to quantum mechanics. Notice that the 
(m, n) — (—1, 0) mode is massless, and that given a localized massless mode, there is also an infinite Kaluza- 
Klein tower of localized massive modes. 

4.3 4D Effective Fermion Charges 

Let us first end this section on fermion fluctuations by obtaining their 4D effective gauge 
couplings. This can be calculated by going beyond their bilinear Lagrangian, and consider- 
ing the interaction term: 

J <fX^G XY M D M X = ... + J d & X^G AT^ + eV^X + ... . (4.40) 

Using the results for the Kaluza-Klein decomposition found in the preceding sections 
X(X) = £ X mn (x)fW(p)e^ = £ X mn (x)^l(u)e- A - B /*e^ , 

m,n m,n 

V»(X) = J2V»rnn(x)f^(py m * = J2 V »rnn(xW^ (4.41) 
m,n m,n 

and recalling that the gauge field zero mode is /qq = 1, a fermion mode X mn has the 
following coupling to the 4D effective gauge group: 



_ e J dipdu tp^nl f$P ip&l 

J dipdu Ipml Ipml 

= e. (4.42) 

Since the gauge field zero mode has a constant wave profile, the effective charges for the 
fermion modes are universal. This is a general result, and independent of any possible 
localization properties of the fermion modes: massless and massive fermion modes will 
always have the same coupling to the massless gauge fields. Again, we find that a large 
mass gap is required in the fermion spectrum in order to hide the Kaluza-Klein tower. We 
now consider this issue in more detail. 
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5 Large Volume Compactifications with a Large Mass Gap 



In the previous two sections, we have calculated the complete Kaluza-Klein spectrum for 
the warped brane world compactification of 6D supergravity, for two interesting sectors of 
the gauge and matter fluctuations. We are now ready to consider the possible implications 
of our results. 

6D brane world models have long been of interest in the context of Large Extra Dimen- 
sions, since these may help with the gauge hierarchy problem. In the conventional ADD 
picture [20], Standard Model particles must be confined to a 4D brane world, in order to 
explain why the large extra dimensions have escaped detection. It would certainly be of 
interest to develop a dynamical description of this localization, within the context of low 
energy effective field theory. 

This could be achieved, for instance, if the zero mode wave profiles were peaked near to a 
brane, and the heavy modes suppressed there [23j. However, we have found that zero mode 
fermions can be peaked near to negative tension branes, only at the price of localizing the 
whole Kaluza-Klein tower (see Figure Therefore, strong couplings are expected between 
light and heavy modes. If the zero mode bulk fermions are to be interpreted as matter in the 
SM, then apparently the only way to explain why we do not observe all the Kaluza-Klein 
modes is by insisting that their mass gap is larger than the WOGeV scale probed to date. 

Usually, this would bring us back to the classical Kaluza-Klein scenario, with the extra 
dimensions required to be very small (at least (lOOGey) -1 scale), and the generation of the 
gauge hierarchy lost. However, in our framework we have seen that a large mass gap can 
occur even if the volume V2, defined as the ratio k 2 / ' k\, becomes large. In the fermionic 
sector this set up can be achieved with the conical defect associated with the warping 
(u) 7^ 1 and ZJ = 1), by taking the parameter r\ to be small. Another way is turning on the 
other defect (ZJ ^ 1) and then taking the large a limit, that is small ZJ limit, corresponding 
to a negative tension brane. The volume V2 in (|2.14|) becomes large but, for both the 
fermions and gauge fields, the mass gap does not reduce to zero. We observe that the latter 
mechanism works also for tq = n, that is uj = ZJ, which corresponds to the unwarped "rugby 
ball" solution with branes of equal tension at each of the two poles |o21 12] . 

The general idea of relaxing the phenomenological constraints on the size of the extra 
dimensions by deforming the shape of the internal space was proposed in |33| . There, it was 
shown that the presence of shape moduli can imply that there is no experimental limit on 
the size of the largest extra dimension. However, requiring a large Kaluza-Klein mass gap 
still constrained the overall volume of the extra dimensions. Here we give an explicit model 
which allows arbitrary large values for both V2 and Mq AP , at least for the fermions and 
vectors. This could have an interesting application in the large extra dimension scenario 
because we can have both k ~ TeV~ 2 and small effects from the massive modes by setting 
a large enough value of Mqap- 

In terms of hiding Large Extra Dimensions from our four dimensional universe, another 
possible approach is to interpret all the bulk fields that we have found (massless modes 
and massive ones) as a hidden sector, only gravitationally coupled to the SM. At this level, 
the SM must then be introduced by hand, confined on the delta-function brane. It seems 
that this is the approach to take if embedding the Supersymmetric Large Extra Dimenions 
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scenario in our calculations, proposed in |2j, to resolve the Cosmological Constant Problem. 
This proposal relates the hierarchy in the Electroweak scale with that of the Cosmological 
Constant. The Electroweak scale is set by the size of the extra dimensions, r, and the 
Cosmological Constant is given by the Kaluza-Klein mass gap, here fixed by the same 
scale 19 1/r. Both may have their observed values when the 6D fundamental scale is TeV, 
and r ~ 0.1mm. For the mass gap to be this small, SM particles must be localized to the 
brane. 

Let us end by considering the tunings involved, when constructing a model with large 
volume (say, y/V% ~ 0.1mm) and large mass gap (say, Mqap ~ TeV^ 1 ). Consider first a 
large mass gap for the fermions. If we set a ~ 1 and r\ <C ro, then the Dirac quantization 
(f2~TT|) implies: 

e !ll ^ jV (5.43) 
ro 91 

If we then assume e ~ 1 (which is natural from group theory) and N ~ 1 (which is required 
for a small number of families) , the large volume - large mass gap condition requires a large 
hierarchy in the bulk gauge couplings: 

~ 10 15 . (5.44) 

9i 

Alternatively, we could set t\ ~ ro and a S> 1, allowing a large mass gap for both 
fermions and gauge fields. In this case, requiring a large mass gap, Mqap ~ V r o> as well 
as large volume, yfV^, ~ \/aro, requires a ~ 10 30 . Then, again, the Dirac quantization 
condition (|2.11|) reveals a large hierarchy in the bulk gauge couplings: 

~ 1(T 30 (5.45) 

91 

In both scenarios we cannot embed the background monopole in U(1)r. 

There are of course other combinations, for example with both a 3> 1 and r\ <C ro, 
in which these hierarchies may be relaxed. However, we should say that these tunings do 
not appear to be very natural or promising. For example, choosing the large dimensionless 
number a 3> 1 corresponds to heavy negative tension branes, and deficit angles orders 
of magnitude less than zero. On the other hand, independently of trying to embed the 
Large Extra Dimension scenario into the present model, we have found an explicit example 
in which the Kaluza-Klein mass gap does not go to zero as the volume goes to infinity, 
contrary to standard lore. 



6 Conclusions 

In this paper we have analyzed an interesting subsector of gauge field and fermion fluc- 
tuations, in the warped brane world solutions of 6D minimal gauged super gravities. In 

19 However, the breakdown of SUSY in the bulk, as in the solutions studied here, may lead to a larger 
prediction for the Cosmological Constant, see 
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Gauge Fields 


Fermions 


V~ (45,l)o 
+(16, l)o 
+(1, 133)o 
+(l,l)o 


A~ (45,1)! 

+(16 + 16, l)i 
+(1,133)! 
~ (l,912)o 



Table 1: The gauge and fermion fields whose Kaluza- Klein spectrum is given by our work, for the illustrative 
example of the anomaly free model Eq x E7 x U(1)r, when the monopole is embedded in E§. We give the 
quantum numbers under Ti = SO(10) x £7 x U(1)r, which is the unbroken subgroup of the 6D gauge group. 

particular, we have focused on bulk components which could give rise to SM or Grand 
Unified gauge and charged matter fields. 

We performed a Fourier decomposition of 6D fields, and transformed the resulting field 
equations into a Schroedinger-like problem. We were then able to find the exact solutions for 
the Kaluza-Klein modes, in terms of hyper geometric functions. We considered in detail the 
boundary conditions that the physical modes must satisfy. In addition to the normalizability 
constraint, consistency also required a hermiticity condition, which can be interpreted as 
demanding current conservation. We were able to implement this in its general, quadratic 
form. Together, these conditions selected the physical modes, and gave rise to a discrete 
mass spectrum, which we presented in full. The discreteness of the spectrum is of course 
to be expected given the compact topology of the internal manifold, whose Euler number 
is two. 

Our study can be applied to several sectors of the 6D supergravities. In Table ^ we 
summarise the 6D fields that are covered by our analysis, for the illustrative example of 
the anomaly free model E§ x Ej x U(l). Moreover, the corresponding spectra for the non- 
super symmetric model of [HI] (at least for the unwarped 4D Poincare invariant case) and 
|32j . generalized to Einstein- Yang- Mills with fermions, can be straightforwardly extracted 
from those given above simply by setting the warp factor to one, that is r\ = r$. 

The exact results presented in this paper enabled us to study the effects of the conical 
defects, sourced by codimension two branes, on the Kaluza-Klein wave profiles and mass 
gaps. As usual, the gauge fields have a zero mode with constant wave profile. For the 
fermions, we found that some zero modes can be peaked on a negative tension brane, but in 
this case the whole Kaluza-Klein tower is peaked there too. Therefore, in order to interpret 
the bulk zero modes as 4D effective fields of the SM, the mass gap must be large. 

Intriguingly, this does not necessarily drive us to the conventional Kaluza-Klein picture, 
with small compact dimensions. It does not, because the conical defects allow a novel 
behaviour in the mass gap, which can be decoupled from the volume of the compactification, 
defined by k 2 / n% . This continues to be observed in the unwarped limit, where the rugby ball 
model of |321 12] is retrieved. Contrary to standard lore, a finite mass gap can be obtained, 
even if the volume goes to infinity. 

For example, a large volume could be arranged in order to generate the Electroweak 
hierarchy, whilst maintaining a large mass gap between the zero modes and the Kaluza-Klein 
tower. This picture does not seem to provide a realisation of the Supersymmetric Large 
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Extra Dimension scenario, where the volume and mass gap should be related. However, in 
this way, SM fields could arise from bulk fields, and along with gravity propagate through 
the large extra dimensions, perfectly consistently with observation. Moreover, for better or 
worse, this picture seems to render the Large Extra Dimension scenario less falsifiable than 
previously thought, since we do not have to expect that the bulk Kaluza-Klein modes are 
accessible at TeV scales. However, arranging for both a large volume and large mass gap 
seems to require a large degree of fine tuning in bulk couplings. Furthermore, for a more 
complete idea, we would have to consider the Kaluza-Klein spectrum for the remaining 
bulk sectors, and in particular the gravitational fluctuations to know the effect of large 
extra dimensions on post-Newtonian tests. 

Indeed, our analysis of 4D effective gauge fields and charged fermions is only a first 
step towards a complete analysis of fluctuations about the warped brane world background 
in 6D supergravity. A final objective would be to derive the full 4D effective field theory 
describing light fluctuations, and an understanding of when 6D physics comes into play. 
The sectors that we have studied here are the simplest ones, in terms of the mixings, but it 
should be possible to continue the project to other fields by extending our work. 

Of the remaining sectors to be analyzed, the scalar perturbations have special impor- 
tance, since they can contain Higgs fields and can have implications for the stability of the 
background solution. Indeed, in the round sphere limit of the model that we have studied, 
tachyons in general emerge from the internal components of the 6D gauge field orthogonal to 
the gauge field background [SSI • A first step in the study of the scalar fluctuations has been 
presented in [12], and it would certainly be interesting if general results can be found. The 
difficulty may be in the mixings of different scalar fluctuations, which lead to a complicated 
system of coupled ODEs. 

In another direction, much of our analysis was general, and could also be used to study 
other theories and other backgrounds with 4D Poincare-2D axial symmetry. 

Finally, it would be interesting to investigate whether there exist other mechanisms, 
which lead to the same decoupling between the mass gap and internal volume that we have 
found here. Indeed, in our set up we have been able to show that the decoupling arises due 
to the conical defects, but it may be possible to find other sources in different frameworks. 
In this way, our explicit example may be a realisation of a more general mechanism. 
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Appendix 



A Conventions and Notation 

We choose the signature —,+,+,+,... for the metric Gmn- The Riemann tensor is defined 
as follows 

R-MNS = ^M^NS ~ ^N^MS + ^MP^NS ~ ^NP^MSi (A-l) 

where the T's are the Levi-Civita connection. Whereas the Ricci tensor and the Ricci scalar 

Rmn = RpMNi R = G MN Rmn- (A. 2) 

Here M,N, ... run over all space-time dimensions. 

Our choice for the 6D constant gamma matrices T A , A = 0, 1, 2, 3, 5, 6, is 

M;:). r*-(y;)> ^{ny ^ 

where the 7^ are the 4D constant gamma matrices and 7 s the 4D chirality matrix. Moreover 
the spin connection is 

< B1 = r, BC Q M c = V B ° (4rM + e A N d M e N B ) , (A.4) 

where ej^ is the vielbein. 

B Delta-Function Singularities 

In this appendix we briefly review how the Ricci scalar acquires a delta-function contribution 
in the presence of a deficit angle, and examine what this implies for our choice of metric 
function e B in (|2.3ft . 

Let us consider the ansatz ()2.3I2.6|) . and illustrate the case for the deficit angle 5 at 
p = 0. Near p = the metric ds 2 , of the 2D internal space can be written as follows 

dsl = d P 2 +(l-^) 2 pV- (B.l) 

By using the change of coordinate r 1-5 / 27r /(l — 5/2ir) = p, this metric becomes 

ds 2 = r- 5 ^ (dr 2 + r 2 dtp 2 ^ . (B.2) 

From (|B.2|) one can show [HE] 

R = 25r s /"5W (y) + (B.3) 
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where the 2D vector y is defined by y = (r cos tp, r sin <p), 5^ is the 2D Dirac delta- function 
and the dots are the smooth contributions 20 . On the other hand, near to p = the Ricci 
scalar, R, can be expressed in terms of derivatives of B: 

R=-B"-\{B'f, (B.4) 

where ' = d p , and from (|B.3|) and ()B.4|) it follows that 

B" = -2 5r s ^8^ (y) + .., (B.5) 

That is, the metric function e B contains a delta-function contribution in its second order 
derivative with respect to p. In coordinate system (|2.4|) . ()B.5|) becomes 

d 2 u B = -2 8r s ^5^ (y) + ..., (B.6) 

The delta-function in the curvature also gives rise to a delta-function in the derivative 
of the spin connection ()4.6|) , The Riemann tensor is defined in terms of the spin connection 
as: 

R A B = dn A B + n A c An c B . (b.7) 

Near to the brane p = 0, relation (|B.7|) gives R p ^ 6 = d p Q^ 6 , and from the expression for 
the spin-connection (|4,6(1 : 

n' = - Qs" + I(£') 2 ^ e B / 2 , (B.8) 

thus leading to the delta- function behaviour from (|B,5|) . 

These results must be recalled when obtaining the Schroedinger-like equations that 
govern the fluctuations. 



C Imposing Boundary Conditions 

Here we study the implications of the NC and the HC for gauge field fluctuations. To 
impose the boundary conditions the following properties will be useful: 

F(a,b,c,z) 1, (C.l) 
F(a, b, c, z) = TiF(a, b, a + b — c + 1, 1 — z) 

+r 2 (l - z) c ~ a ~ h F{c - a,c - b, c- a - b + 1, 1 - z), (C.2) 

where 

_ r(c)r(c-o-6) _ T(c)T(-c + a + b) 

_^_^c_-a)T{c-b) , 2 ~ r(a)T(6) ' { ' 

20 Eq. <B.3t describes the asymptotic behaviour in the vicinity of the brane |36|. It seems that for positive 
S the Ricci scalar vanishes at the origin. However, the first term on the right hand side of (B.3) should be 
interpreted in a distributional sense. The effect of the Ricci scalar as a distribution on a scalar test function 
/(y) is then J d 2 y^gRf(y) = 2 5/(0) + . . .. 
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and r is the Euler gamma function. The relation (jC.2|) is valid if c — a — b is not an integer 
37 and y\ and y 2 in (|3,21|) are both well defined when c is not an integer. In general c—a — b 
and c are not integers for generic u and U; so we can consider u and uJ as regulators to use 
(|3.21|) and (|C2|) and at the end we can take the limits in which c — a — b and c go to an 
integer, which will turn out to be well defined. 

We first consider the behaviour of tp for u — ► u, that is z — ► because of the definition 
z = cos 2 (j^j- F° r c / 1 we use the expression for ift given in (|3.22j) and property (|C.1|) 
gives us 

ij; u =T Kxiu-u) 2 ^ + K 2 {u-u) 1 - 2 ^, (C.4) 

where we used c = 1/2 + 27. So the NC (|3.5|) implies K 1 = when 7 < —1/4 and K2 = 
when 7 > 3/4. On the other hand the HC ()3.16|) implies 21 

lim V* (~d u + l^-z) rp < 00 (C.5) 
lt-»u \ 2 u — u/ 

and by using the behaviour ()C.4|) this limit becomes 

(2^ - I) lim [|^i| 2 (n - n) 4 ^ 1 - K\K 2 + K X K* 2 - \K 2 \ 2 (u - u)-^ +1 ] , (C.6) 

\ 2 J u-^u L J 

so the HC implies K\ = when 7 < 1/4 and K 2 = when 7 > 1/4. The case 7 = 1/4 
corresponds to c = 1 and so we have to use the expression of tp given in Q3.24JI . We have 
then 

V> u =? Kiiu - ufl 2 -K 2 (u- u) 1 / 2 ln(n - u) (C.7) 

for which (|C.5|) implies K 2 = 0. Therefore we obtain (|3.25j) and (|3,26|) , 

The discreteness of the spectrum emerges when we impose the NC and HC for u — > 0. 
For instance for m > Ny, up to an overall constant, the behaviour of ip is given by properties 
jSIH and (iri2l : 

^4°^^ + ^, (C.8) 

where 2 are defined in (|C.3|) and we used c — a — b = 1/2 — 2/3. Behaviour (|C.8|) is similar 
to (|C"4ll but 7 is replaced by /3. So, following the same steps as above, the NC and the HC 
imply that T\ = for /3 < 1/4 and T 2 = for /5 > 1/4. Let us study the case m > Ny and 
< 1/4, that is 

AV < m < 0. (C.9) 



We then have 



T(c — a)T(c — b) 



21 In fact, for the gauge field sector, each term in the HC 13.161 is separately zero if one requires their 
finiteness. Therefore l|C,50 and its counter-part for u = are sufficient to ensure l|3.16|l . On the other hand, 
for the fermions, there are some cases in which they are each finite and non-zero, and so l|3.16t requires that 
they cancel. 
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Since the Euler gamma function never vanishes we require that T(c— a) = oo or T(c—b) = oo 
and this is possible only when c — a = —n or c — b = —n, where n = 0, 1, 2, 3, .... By using 
the definitions (|3.19j) both conditions lead to the following squared masses 



M 2 

V n,m 



2 |n(n + 1) + Q + rij \-mu) + (m - iVy)aJ] J (C.ll) 

which are positive because of ()C.9|) and n > 0. When m > iVy and 22 /3 > 1/4, that is 

m > N v and m> 0, (C.12) 

we have 

_ T(c)T(-c + a + b) 

°- T2 = — fww) — ( } 

and this implies a = —n or b = —n, where n = 0,1,2,3,.... The corresponding squared 
masses are 



M 2 

1V1 V n,m 



\ |n(n + 1) + Q +nj [mco + (m - iV»aJ] + mw(m - iVy)uj| (C.14) 



which are positive or vanishing. We can study the case m < Ny in a similar way. The 
complete result for the gauge fields sector is given in equations ()3.27j) - ()3.3fl|) . We observe 
that (|3.16|) is now automatically satisfied by every pair of wave functions ip and ip', for a given 
quantum number m, since the asymptotic behaviour of the wave function cannot depend 
on the quantum number n: this is a consequence of the 1/it 2 and l/(u — u) 2 singularities 
of the potential V in (|3.14f) . 



D Complete Fermionic Mass Spectrum 

In this appendix we give the complete fermionic spectrum which is also labeled by an integer 
quantum number n = 0, 1, 2, .... Although much longer, the calculation proceeds in exactly 
the same way as for the gauge field sector, outlined in the previous appendix 23 . 



m>~l + N + ^ 



in this case = and we obtain the following squared masses. 



For m > -\ + ± 



M 



F n,m 





n v 




^ + u + lv 


— h m 






V2 J. 





1 -f 1 

— hn + w — \-m 

2 V2 



N 



> 0. 



(D.l) 



22 The (3 = 1/4 case is recovered by taking the limit u> — *• 0. 

23 There is one additional subtlety. Here, for the values of m which allow a zero mode, we must impose a 
mixed HC between the massless mode and massive modes, in addition to the diagonal HC. In general the 
HC involving distinct wave functions ipmn and tp mn i does not lead to additional constraints, because the 
asymptotic behaviour of the modes is independent of n. However, the massless modes are more strongly 
constrained than the massive ones, obeying as they do a decoupled Dirac equation in addition to the 
Schroedinger equation. 
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M 



Fn,m 



1 



- + n + lo — + m 



1 



- + n + uj[ - + m-N 



> 0. (D.2) 



or 



Mi 



Fn,m 



r(l + n) 



1 



1 + n + w - + m — AM — w -+m 



1 



> 0. (D.3) 



• For m < -i - 



M 



Fn,m 



1 



l + n + cj - + m — A 7 ] — a; - + m 



1 



> 0. (D.4) 



?7l<— ^ + A^— 2= : in this case K\ = and we obtain the following squared masses. 



Form>-i + 2b 



r 



n — oj [ — h m — TV J + lo — \-m 



> 0. 



(D.5) 



• Vox -\-±<m<-\ + ± 



M Fn, m = ~2 n 

r 



n — to [ - + m — N ) + lo ( - + m 



> 0. 



(D.6) 



• Form< 



M 



Fn,m 



+ n — LO \ — + m 



+ n — uj\- + m — N 



> 0. (D.7) 



— ^ + A^— 2=<?7l<— ^ + A^+2= : this case is possible only when 5 < 0. 
• For m > — \ + we have Ki = and 



M 



Fn,m 



Til 



n — lo { - + m — N ) + lo ( - + m 



> 0. 



(D.8) 



For m < — \ — we have two possibilities. We have Ki = and 



M ln, m = 12(1 + ^) 



1 + n + cj — \-m — N } — to I — h m 



1 



> (D.9) 



or K\ = and 



M 2 - — 

M Fn,m - 2 
r 



+ n — lo \ — Yra 



+ n — lo \ — \- m — N 



> (D.10) 
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_ 1 

2 

#1 



2w 

and 



< m < 



.1 + _L. 

2 ^ 2uT 



this case is possible only when 5 < and we obtain 



M 



F n,m 



n 



+ m — N 



+ <jj \ - + m 



> 0. 



(D.ll) 



Again, we can perform a check of our results by considering the S 2 limit (uj — ► 1, u — ► 1). 
In this case, the mass spectrum ()D.1|) - ()D.11|) reduces correctly to 

a 2 Ml =(l + 1 + — ) (l + 1 - — ) , multiplicity = 21 + 1 (D.12) 



where 24 / = ^1 1 + A; and k = 0, 1, 2, 3, .... 
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